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Abstract 

An algorithm M is described that solves any well-defined problem p as quickly 
as the fastest algorithm computing a solution to p, save for a factor of 5 and low- 
order additive terms. M optimally distributes resources between the execution 
of provably correct p-solving programs and an enumeration of all proofs, including 
relevant proofs of program correctness and of time bounds on program runtimes. M 
avoids Blum's speed-up theorem by ignoring programs without correctness proof. M 
has broader applicability and can be faster than Levin's universal search, the fastest 
method for inverting functions save for a large multiplicative constant. An extension 
of Kolmogorov complexity and two novel natural measures of function complexity 
are used to show that the most efficient program computing some function / is also 
among the shortest programs provably computing /. 



1 Publishcd in the International Journal of Foundations of Computer Science, Vol. 13, No. 3, (June 
2002) 431-443. Extended version of An effective Procedure for Speeding up Algorithms (cs.CC/0102018) 
presented at Workshops MaBiC-2001 & TAI-2001. 
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1 Introduction & Main Result 

Searching for fast algorithms to solve certain problems is a central and difficult task in 
computer science. Positive results usually come from explicit constructions of efficient 
algorithms for specific problem classes. A wide class of problems can be phrased in the 
following way. Given a formal specification of a problem depending on some parameter 
x G X, we are interested in a fast algorithm computing solution y 6 Y. This means 
that we are interested in a fast algorithm computing / : X — > Y, where / is a formal 
(logical, mathematical, not necessarily algorithmic), specification of the problem. Ideally, 
we would like to have the fastest algorithm, maybe apart from some small constant factor 
in computation time. Unfortunately, Blum's Speed-up Theorem ||, Q shows that there 
are problems for which an (incomputable) sequence of speed-improving algorithms (of 
increasing size) exists, but no fastest algorithm. 

In the approach presented here, we consider only those algorithms which provably solve 
a given problem, and have a fast (i.e. quickly computable) time bound. Neither the 
programs themselves, nor the proofs need to be known in advance. Under these constraints 
we construct the asymptotically fastest algorithm save a factor of 5 that solves any well- 
defined problem /. 

Theorem 1. Let p* be a given algorithm computing p* (x) from x, or, more generally, a 
specification of a function. Let p be any algorithm, computing provably the same function 
as p* with computation time provably bounded by the function t p (x) for all x. time tp (x) 
is the time needed to compute the time bound t p (x). Then the algorithm M p * constructed 
in Section [| computes p*(x) in time 

timeM p * (%) < b-t p (x) + d p -time tp (x) + c p 

with constants c p and d p depending on p but not on x. Neither p, t p , nor the proofs need 
to be known in advance for the construction of M p * (x) . 

Known time bounds for practical problems can often be computed quickly, i.e. 
timet p (x)/time p (x) often converges very quickly to zero. Furthermore, from a practi- 
cal point of view, the provability restrictions are often rather weak. Hence, we have 
constructed for every problem a solution, which is asymptotically only a factor 5 slower 
than the (provably) fastest algorithm! There is no large multiplicative factor and the 
problems are not restricted to inversion problems, as in Levin's algorithm (see section ||]). 
What somewhat spoils the practical applicability of M p * is the large additive constant c p , 
which will be estimated in Section [|. 

An interesting and counter-intuitive consequence of Theorem 1, derived in Section |7|, is 
that the fastest program that computes a certain function is also among the shortest 
programs that provably computes this function. Looking for larger programs saves at 
most a finite number of computation steps, but cannot improve the time order. 
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In section ^| we review Levin search and the universal search algorithms simple and 
search, described in [13|]. We point out that simple has the same asymptotic time 



complexity as SEARCH not only w.r.t. the problem instance, but also w.r.t. to the prob- 
lem class. In Section |3] we elucidate Theorem 1 and the range of applicability on an 
example problem unsolvable by Levin search. In Section |] we give formal definitions of 
the expressions time, proof, compute, etc., which occur in Theorem 1, and define the fast 
algorithm M p * . In Section [5] we analyze the algorithm M p * , especially its computation 
time, prove Theorem 1, and give upper bounds for the constants c p and d p . Subtleties 
regarding the underlying machine model are briefly discussed in Section || In Section [7] 
we show that the fastest program computing a certain function is also among the shortest 
programs provably computing this function. For this purpose, we extend the definition 
of the Kolmogorov complexity of a string and define two new natural measures for the 
complexity of functions and programs. Section § outlines generalizations of Theorem 1 to 
i/o streams and other time- measures. Conclusions are given in Section BL 



2 Levin Search 

Levin search is one of the few rather general speed-up algorithms. Within a (typically 
large) factor, it is the fastest algorithm for inverting a function g : Y — > X, if g can 
be evaluated quickly (JTTJ, 0. Given x, an inversion algorithm p tries to find a y £ Y, 



called g-witness for x, with g(y) =x. Levin search just runs and verifies the result of all 
algorithms p in parallel with relative computation time 2~ 1 ^; i.e. a time fraction 2~ 1 ^ is 
devoted to execute p, where l(p) is the length of program p (coded in binary). Verification 
is necessary since the output of any program can be anything. This is the reason why 
Levin search is only effective if a fast implementation of g is available. Levin search halts 
if the first g-witness has been produced and verified. The total computation time to find a 
solution (if one exists) is bounded by 2 l( - p ^-time p (x). time p (x) is the runtime of p(x) plus 
the time to verify the correctness of the result {g{p{x)) =x) by a known implementation 
for g. 

Li and Vitanyi Jl3|, p503] propose a very simple variant, called simple, which runs all 
programs P1P2P3 ■ ■ ■ one step at a time according to the following scheme: pi is run every 
second step, p 2 every second step in the remaining unused steps, p 3 every second step 
in the remaining unused steps, and so forth, i.e. according to the sequence of indices 
121312141213121512 .... If pk inverts g on x in time Pk (x) steps, then simple will do the 
same in at most 2 k time pk (x) + 2 fe_1 steps. In order to improve the factor 2 k , they define 
the algorithm search, which runs all p (of length less than i) for 2 l 2~ l ^ p > steps in phase 
% = 1, 2, 3, . . ., until it has inverted g on x. The computation time of search is bounded 
by 2 K( - k ^ + °^ l Hime pk {x), where K{k) <l{pk) <21og/c is the Kolmogorov complexity of k. 
They suggest that simple has worse asymptotic behaviour w.r.t. k than search, but 
actually this is not the case. 

In fact, simple and search have the same asymptotics also in k, because search itself 
is an algorithm with some index fc SEARCH =0(1). Hence, simple executes search every 
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2 fcsEARCH -th step, and can at most be a constant (in k and x) factor 2 fcsEARCH = 0(1) slower 
than search. However, in practice, search should be favored, because also constants 
matter, and 2 fcsEARCH is rather large. 



Levin search can be modified to handle time-limited optimization problems as well j2C 
Many, but not all problems, are of inversion or optimization type. The matrix multi- 
plication example (section [^), the decision problem SAT [ 13], p503], and reinforcement 



learning ||, for instance, cannot be brought into this form. Furthermore, the large factor 
2'(p) somewhat limits the applicability of Levin search. See |T3| , pp518-519] for a historical 
review and further references. 

Levin search in program space cannot be used directly in M p * for computing p* because 
we have to decide somehow whether a certain program solves our problem or computes 
something else. For this, we have to search through the space of proofs. In order to 
avoid the large time-factor 2 l ( p \ we also have to search through the space of time-bounds. 
Only one (fast) program should be executed for a significant time interval. The algo- 
rithm M p * essentially consists of 3 interwoven algorithms: sequential program execution, 
sequential search through proof space, and Levin search through time-bound space. A 
tricky scheduling prevents performance degradation from computing slow p's before the p 
has been found. 



3 Applicability of the Fast Algorithm M p * 

To illustrate Theorem 1, we consider the problem of multiplying two n x n matrices. If 
p* is the standard algorithm for multiplying two matrices^ xER nn xR n " n of size l(x)~n 2 , 
then t p *(x) :=2n 3 upper bounds the true computation time time p *(x) = n 2 (2n — 1). We 
know there exists an algorithm p' for matrix multiplication with time p i(x) < t p /(x) : = 
c-n 2 ' 81 [21|. The time-bound function (cast to an integer) can, as in many cases, be 
computed very quickly, time t ,{x) = 0{log 2 n). Hence, using Theorem 1, also M p * is fast, 
timeM p * (%) < 5en 2 ' 81 + 0(log 2 n) . Of course, M p * would be of no real use if p' is already the 
fastest program, since p' is known and could be used directly. We do not know however, 
whether there is an algorithm p" with time p n(x) < d-n 2 logn, for instance. But if it does 
exist, timeM p *(x) < 5d-n 2 log n+0(l) for all x is guaranteed. 

The matrix multiplication example has been chosen for specific reasons. First, it is not 
an inversion or optimization problem suitable for Levin search. The computation time 
of Levin search is lower-bounded by the time to verify the solution (which is at least 
c-ra 2 ' 81 to our knowledge) multiplied with the (large) number of necessary verifications. 
Second, although matrix multiplication is a very important and time-consuming issue, p' 
is not used in practice, since c is so large that for all practically occurring n, the cubic 
algorithm is faster. The same is true for c p and d p , but we must admit that although c 
is large, the bounds we obtain for c p and d p are tremendous. On the other hand, even 
Levin search, which has a tremendous multiplicative factor, can be successfully applied 

2 Instead of interpreting R as the set of real numbers one might take the field F2 ~ {0, 1} to avoid 
subtleties arising from large numbers. Arithmetic operations are assumed to need one unit of time. 
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14] , |i~6|| , when handled with care. The same should hold for Theorem 1, as will be 
discussed. We avoid the OQ notation as far as possible, as it can be severely misleading 
(e.g. 10 42 = 0(1)°^ = 0(1)). This work could be viewed as another OQ warning 
showing, how important factors, and even subdominant additive terms, are. 

An obvious time bound for p is the actual computation time itself. An obvious algorithm 
to compute time p (x) is to count the number of steps needed for computing p(x). Hence, 
inserting t p = time p into Theorem 1 and using time tirnep (x) <time p (x), we see that the 
computation time of M p * is optimal within a multiplicative constant (d p + 5) and an 
additive constant c p . The result is weaker than the one in Theorem 1, but no assumption 
concerning the computability of time bounds has to be made. 

When do we trust that a fast algorithm solves a given problem? At least for well specified 
problems, like satisfiability, solving a combinatoric puzzle, computing the digits of tt, 
we usually invent algorithms, prove that they solve the problem and in many cases also 
can prove good and fast time bounds. In these cases, the provability assumptions in 
Theorem 1 are no real restriction. The same holds for approximate algorithms which 
guarantee a precision e within a known time bound (many numerical algorithms are of 
this kind). For exact /approximate programs provably computing/converging to the right 
answer (e.g. traveling salesman problem, and also many numerical programs), but for 
which no good, and easy to compute time bound exists, M pjf is only optimal apart from 
a huge constant factor 5 + d p in time, as discussed above. A precursor of algorithm M p * 
for this case, in a special setting, can be found in ||Q. For poorly specified problems, 
Theorem 1 does not help at all. 



4 The Fast Algorithm M p * 

One ingredient of algorithm M p * is an enumeration of proofs of increasing length in some 
formal axiomatic system. If a proof actually proves that p and p* are functionally equiv- 
alent and p has time bound t p , add (p,t p ) to a list L. The program p in L with the 
currently smallest time bound t p (x) is executed. By construction, the result p(x) is iden- 
tical to p*(x). The trick to achieve the time bound stated in Theorem 1 is to schedule 
everything in a proper way, in order not to lose too much performance by computing slow 
p's and t p s before the p has been found. 

To avoid confusion, we formally define p and t p to be binary strings. That is, p is neither 
a program nor a function, but can be informally interpreted as such. A formal definition 
of the interpretations of p is given below. We say "p computes function f", when a 
universal reference Turing machine U on input (p, x) computes f(x) for all x. This is 
denoted by U (p, x) = f(x). To be able to talk about proofs, we need a formal logic system 
(V, A, i/i, Ci, fi, Ri, — >, A, =, ...), and axioms, and inference rules. A proof is a sequence of 
formulas, where each formula is either an axiom or inferred from previous formulas in the 

3 The algorithm AI£* J creates an incremental policy for an agent in an unknown non-Markovian 
environment, which is superior to any other time t and space I bounded agent. The computation time of 
AI£ tz is of the order t-2 l . 
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sequence by applying the inference rules. See || or any other textbook on logic or proof 
theory. We only need to know that provability, Turing Machines, and computation time 
can be formalized: 

1. The set of (correct) proofs is enumerable. 

2. A term u can be defined such that the formula \iy:u{p,y)=u{p*,y)} is true if, and 
only if U(p, x) = U(p*, x) for all x, i.e. if p and p* describe the same function. 

3. A term tm can be defined such that the formula [tm(p, x) = n] is true if, and only if 
the computation time of U on (p, x) is n, i.e. if n = time p (x). 

We say that p is provably equivalent to p* if the formula [Wy : u(p, y) = u{p* , y)\ can be 
proved. M p * runs three algorithms A, B, and C in parallel: 

Algorithm M p *(x) 

Initialize the shared variables L : = {}, tf ast := oo, Pf as t '■= V* ■ 
Start algorithms A, B, and C in parallel with 10%, 10% and 80% 
computational resources, respectively. 

That is, C performs 8 steps when A and B perform 1 step each. 

Algorithm A 

for i:=l,2,3,... do 

pick the i th proof in the list of all proofs and 
isolate the last formula in the proof. 

if this formula is equal to \\/y:u(p, y) = u(p*, y) A u(t, y) > tm{p, y)} 
for some strings p and t, 
then add {p, t) to L. 
next i 

Algorithm B 

for all (p,t)eL 

run U on all (t, x) in parallel for all t with relative computational resources 2~ l ( p ^~ l( - t \ 
if U halts for some t and U(t,x) <tf ast , 
then t fast ■= U(t,x) and p fast := p. 
continue (p, t) 

Algorithm C 

for k:=l,2,4,8,16,32,... do 

pick the currently fastest program p := Pf ast with time bound tf ast . 
run U on (p, x) for k steps, 
if U halts in less than k steps, 

then print result U(p, x) and abort computation of A, B and C. 
continue k. 

Note that A and B only terminate when aborted by C. Figure [1] illustrates the time- 
scheduling on a fictitious example. The discussion of the algorithm(s) in the following 
sections clarifies details and proves Theorem 1. 
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Figure 1: Sample execution of M p *. The horizontal axis marks the progress of M p *. 
The vertical axis marks general time-like quantities. The lower solid stepped curve is the 
content of variable tf ast at time t. Black dots on this curve (denoted by p) mark events 
where algorithm A creates and adds another valid program p together with its time-bound 
t to list L. The (later occurring) step in the curve marks the time when algorithm B 
has computed t(x) and updates tf as t- Algorithm C starts p at the next instruction step k 
which is a power of 2 and executes p for k steps. This is illustrated by the exponentially 
increasing sawtooth curve. The time bound of p during execution is indicated by the 
upper dashed stepped curve. Algorithm M p * terminates when the stepped curve crosses the 
sawtooth curve (white circle), or earlier if the time-bound is not sharp. 

5 Time Analysis 

Henceforth we return to the convenient abbreviations p(x) :=U(p, x) and t p (x) :—U(t p , x). 
Let p' be some fixed algorithm that is provably equivalent to p*, with computation time 
time p i provably bounded by t p >. Let l(proof(p')) be the length of the binary coding 
of the, for instance, shortest proof. Computation time always refers to true overall 
computation time, whereas computation steps refer to instruction steps, steps = a -time, 
if a percentage a of computation time is assigned to an algorithm. 

A) To write down (not to invent!) a proof requires 0(l(proof)) steps. A time 
0(N axiom - l(Fi)) is needed to check whether a formula Fi in the proof FiF 2 ...F n is an 
axiom, where N axiom is the number of axioms or axiom-schemes, which is finite. Variable 
substitution (binding) can be performed in linear time. For a suitable set of axioms, the 
only necessary inference rule is modus ponens. If Fi is not an axiom, one searches for 
a formula Fj, j < i of the form Fk — > Fi and then for the formula Fk, k <i. This takes 
time 0(l(proof)). There are n < 0(l(proof)) formulas Fi to check in this way. Whether 
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the sequence of formulas constitutes a valid proof can, hence, be checked in 0(1 (proof) 2 ) 
steps. There are less than 2 l+1 proofs of (binary) length < I. Algorithm A receives 
a = 10% of relative computation time. Hence, for a proof of (p',t p >) to occur, and for 
0', V) to be added to L, at most time T A < -^■2 l ^ roo ^ +x -0(l(proof(p')) 2 ) is needed. 
Note that the same program p can and will be accompanied by different time bounds t p ; 
for instance (p, time p ) will occur. 



B) The time assignment of algorithm B to the t p s only works if the Kraft inequality 

i(j>)-i(t p ) 



E 



(p,t p )£L ■ 



< 1 is satisfied [HJ. This can be ensured by using prefix free 

,,(x) is, by 



e.g. Shannon- Fano) codes |L7|, |13| . The number of steps to calculate t 1 



definition, timet , (x) . The relative computation time a available for computing t p /(x) 
is 10% • 2~^')~'(V). Hence, t p i(x) is computed and tf ast < t p i(x) is checked after time 
T B < T A + j5%-2 l(p ' )+l{tp,) ■time tp ,(x). We have to add T A , since B has to wait, in the 
worst case, time T A before it can start executing t p i{x). 



C) If algorithm C halts, its construction guarantees that the output is correct. In the 
following, we show that C always halts, and give a bound for the computation time. 

i) Assume that algorithm C stops before B performed the check t p >(x) < tf ast , because 
a different p already computed p(x). In this case Tc < Tb- 

ii) Assume that k = ko in C when B performs the check t p /(x) < tf ast . Running-time 
Tb has passed until this point, hence ko < 80% -Tb ■ Furthermore, assume that C 
halts in period ko because the program (different from p') executed in this period 
computes the result. In this case, Tc < -^^2ko < 2T B - 

in) If C does not halt in period k but 2k > tf ast , then p'{x) has enough time to compute 
the solution in the next period k = 2ko, since time p >(x) < tf ast < Ako — 2ko. Hence 
Tc < sk^ko < AT B . 

iv) Finally, if 2k < tf ast we "wait" for the period k> k with \k < tf ast < k. In this 
period k, either p'(x), or an even faster algorithm, which has in the meantime been 
constructed by A and B, will be computed. In any case, the 2k — k > tf ast steps are 
sufficient to compute the answer. We have 80% -T c < 2k < 4tf ast < 4t p i(x). 



The maximum of the cases (i) to (iv) bounds the computation time of C and, hence, of 
M r by 

timeM p *(x) = Tc < max{4Te, 5t p (x)} < ATb + 5t p (x) < 

< 5-t p (x) + d p -time tp (x) + c p 

d p = 40-2' (p)+i( H c p = 40-2 l{proof{p))+1 -O(l(proof(p) 2 ) 

where we have dropped the prime from p. We have also suppressed the dependency of 
c p and dp on p* (proof (p) depends on p* too), since we considered p* to be a fixed given 
algorithm. The factor of 5 may be reduced to 4 + e by assigning a larger fraction of time 
to algorithm C. The constants c p and d p will then be proportional to -. We were not 
able to further reduce this factor. 
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6 Assumptions on the Machine Model 

In the time analysis above we have assumed that program simulation with abort possibil- 
ity and scheduling parallel algorithms can be performed in real-time, i.e. without loss of 
performance. Parallel computation can be avoided by sequentially performing all opera- 
tions for a limited time and then restarting all computations in a next cycle with double 
the time and so on. This will increase the computation time of A and B (but not of C!) 
by, at most, a factor of 4. Note that we use the same universal Turing machine U with 
the same underlying Turing machine model (number of heads, symbols, ...) for measur- 
ing computation time of all programs (strings) p, including M p *. This prevents us from 
applying the linear speedup theorem (which is cheating somewhat anyway), but allows 
the possibility of designing a U which allows real-time simulation with abort possibility. 
Small additive "patching" constants can be absorbed in the 0() notation of c p . Theorem 
1 should also hold for Kolmogorov-Uspenskii and Pointer machines. 



7 Algorithmic Complexity and the Shortest Algo- 
rithm 

Data compression is a very important issue in computer science. Saving space or channel 
capacity are obvious applications. A less obvious (but not far fetched) application is that 
of inductive inference in various forms (hypothesis testing, forecasting, classification, ...). 
A free interpretation of Occam's razor is that the shortest theory consistent with past 
data is the most likely to be correct. This has been put into a rigorous scheme by 



and proved to be optimal in jjl9j Kolmogorov Complexity is a universal notion of 



the information content of a string || |, [23| . It is defined as the length of the shortest 
program computing string x. 

K v (x) := mm{l(p) : U(p) = x} = K{x) + 0(1) 

where U is some universal Turing Machine. It can be shown that Kjj(x) varies, at most, by 
an additive constant independent of x by varying the machine U. Hence, the Kolmogorov 
Complexity K{x) is universal in the sense that it is uniquely defined up to an additive 
constant. K(x) can be approximated from above (is co-enumerable), but not finitely 



computable. See [ 13 1 for an excellent introduction to Kolmogorov Complexity and 



for a review of Kolmogorov inspired prediction schemes. 



Recently, Schmidhuber [To] has generalized Kolmogorov complexity in various ways to 
the limits of computability and beyond. In the following, we also need a generalization, 
but of a different kind. We need a short description of a function, rather than a string. 
The following definition of the complexity of a function / 

K'(f) := min{/(p) : U(p,x) = f(x) Vx} 
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seems natural, but suffers from not even being approximable. There exists no algorithm 
converging to K'(f), because it is undecidable whether a program p is the shortest pro- 
gram equivalent to a function /. Even if we have a program p* computing /, K'(p*) is 
not approximable. Using K(p*) is not a suitable alternative, since K(p*) might be con- 
siderably longer than K'(p*), as in the former case all information contained in p* will 
be kept - even that which is functionally irrelevant (e.g. dead code). An alternative is to 
restrict ourselves to provably equivalent programs. The length of the shortest one is 

K"(p*) : = min{/(p) : a proof of \Vy.u(p,y) = u(p*,y)) exists} 

It can be approximated from above, since the set of all programs provably equivalent to 
p* is enumerable. 

Having obtained, after some time, a very short description p' of p* for some purpose 
(e.g. for defining a prior probability for some inductive inference scheme), it is usually 
also necessary to obtain values for some arguments. We are now concerned with the 
computation time of p' . Could we get slower and slower algorithms by compressing p* 
more and more? Interestingly this is not the case. Inventing complex (long) programs 
is not necessary to construct asymptotically fast algorithms, under the stated provability 
assumptions, in contrast to Blum's Theorem 0, |3|. The following theorem roughly says 
that there is a single program, which is the fastest and the shortest program. 

Theorem 2. Let p* be a given algorithm or formal specification of a function. There 
exists a program p, equivalent to p* , for which the following holds 

i) l(p) <K"(p*) + 0(l) 
ii) time p (x) < 5-t p (x) + d p -time t (x) + c p 

where p is any program provably equivalent to p* with computation time provably less 
than t p (x). The constants c p and d p depend on p but not on x. 

To prove the theorem, we just insert the shortest algorithm p' provably equivalent to p* 
into M, that is p := M p >. As only 0(1) instructions are needed to build M p i from p', M p i 
has size l(p')+0(l) = K"(p*)+0(1). The computation time of M p > is the same as of M p * 
apart from "slightly" different constants. 

The following subtlety has been pointed out by Peter van Emde Boas. Neither M p * , nor p 
is provably equivalent to p* . The construction of M p * in section |] shows equivalence of M p * 
(and of p) to p*, but it is a meta-proof which cannot be formalized within the considered 
proof system. A formal proof of the correctness of M p * would prove the consistency of 
the proof system, which is impossible by Godels second incompleteness theorem. See 
for details in a related context. 
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8 Generalizations 

If p* has to be evaluated repeatedly, algorithm A can be modified to remember its current 
state and continue operation for the next input (A is independent of x\). The large offset 
time c p is only needed on the first invocation. 

M p * can be modified to handle i/o streams, definable by a Turing machine with monotone 
input and output tapes (and bidirectional working tapes) receiving an input stream and 
producing an output stream. The currently read prefix of the input stream is x. time p (x) 
is the time used for reading x. M p * caches the input and output streams, so that algorithm 
C can repeatedly read/write the streams for each new p. The true input/output tapes are 
used for request ing/producing a new symbol . Algorithm B is reset after 1, 2, 4, 8, ... steps 
(not after reading the next symbol of x\) to appropriately take into account increased 
prefixes x. Algorithm A just continues. The bound of Theorem 1 holds for this case too, 
with slightly increased d p . 

The construction above also works if time is measured in terms of the current output rather 
than the current input x. This measure is, for example, used for the time-complexity of 
calculating the n th digit of a computable real (e.g. 7r), where there is no input, but only 
an output stream. 



9 Summary & Outlook 

We presented an algorithm M p * which accelerates the computation of a program p*. M p * 
combines (A) sequential search through proof space, (B) Levin search through time- 
bound space, (C) and sequential program execution, using a somewhat tricky scheduling. 
Under certain provability constraints, M p * is the asymptotically fastest algorithm for 
computing p* apart from a factor 5 in computation time. Blum's Theorem shows that 
the provability constraints are essential. We have shown that the conditions on Theorem 
1 are often, but not always, satisfied for practical problems. For complex approximation 
problems, for instance, where no good and fast time bound exists, M p * is still optimal, 
but in this case, only apart from a large multiplicative factor. We briefly outlined how 
M p * can be modified to handle i/o streams and other time- measures. An interesting and 
counter-intuitive consequence of Theorem 1 was that the fastest program computing a 
certain function is also among the shortest programs provably computing this function. 
Looking for larger programs saves at most a finite number of computation steps, but 
cannot improve the time order. To quantify this statement, we extended the definition 
of Kolmogorov complexity and defined two novel natural measures for the complexity of 
a function. The large constants c p and d p seem to spoil a direct implementation of M p *. 
On the other hand, Levin search has been successfully applied to solve rather difficult 



machine learning problems [|14], |16[ , even though it suffers from a large multiplicative 
factor of similar origin. The use of more elaborate theorem-provers, rather than brute 
force enumeration of all proofs, could lead to smaller constants and bring M* closer to 
practical applications, possibly restricted to subclasses of problems. A more fascinating 
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(and more speculative) way may be the utilization of so called transparent or holographic 
proofs [|TJ. Under certain circumstances they allow an exponential speed up for checking 
proofs. This would reduce the constants c p and d p to their logarithm, which is a small 
value. I would like to conclude with a general question. Will the ultimate search for 
asymptotically fastest programs typically lead to fast or slow programs for arguments of 
practical size? Levin search, matrix multiplication and the algorithm M p * seem to support 
the latter, but this might be due to our inability to do better. 



Acknowledgements 

Thanks to Monaldo Mastrolilli, Jurgen Schmidhuber, and Peter van Emde Boas for en- 
lightening discussions and for useful comments. This work was supported by SNF grant 
2000-61847.00. 



References 

L. Babai, L. Fortnow, L. A. Levin, and M. Szegedy. Checking computations in 
polylogarithmic time. STOC: 23rd ACM Symp. on Theory of Computation, 23:21- 
31, 1991. 

M. Blum. A machine-independent theory of the complexity of recursive functions. 
Journal of the ACM, 14(2):322-336, 1967. 

M. Blum. On effective procedures for speeding up algorithms. Journal of the ACM, 
18(2):290-305, 1971. 

G. J. Chaitin. On the length of programs for computing finite binary sequences. 
Journal of the ACM, 13(4):547-569, 1966. 

Melvin C. Fitting. First-Order Logic and Automated Theorem Proving. Graduate 
Texts in Computer Science. Springer- Verlag, Berlin, 2nd edition, 1996. 

J. Hartmanis. Relations between diagonalization, proof systems, and complexity 
gaps. Theoretical Computer Science, 8(2):239-253, April 1979. 

M. Hutter. New error bounds for Solomonoff prediction. Journal of Computer and 
System Science, in press, 1999. |ftp: / / ftp.idsia.ch / pub/ techrep/IDSIA-1 l-00.ps.gz[ 



M. Hutter. A theory of universal artificial intelligence based on algorithmic complex- 
ity. Technical report, 62 pages, 2000. |ittp: / / arxiv.org/ abs/ cs. Al/ 0004001 . 



A. N. Kolmogorov. Three approaches to the quantitative definition of information. 

Problems of Information and Transmission, l(l):l-7, 1965. 

L. G. Kraft. A device for quantizing, grouping and coding amplitude modified pulses. 
Master's thesis, Cambridge, MA, 1949. 



Marcus Hutter, The Fastest and Shortest Algorithm 



12 



L. A. Levin. Universal sequential search problems. Problems of Information Trans- 
mission, 9:265-266, 1973. 

L. A. Levin. Randomness conservation inequalities: Information and independence 
in mathematical theories. Information and Control, 61:15-37, 1984. 

M. Li and P. M. B. Vitanyi. An introduction to Kolmogorov complexity and its 
applications. Springer, 2nd edition, 1997. 

J. Schmidhuber. Discovering neural nets with low Kolmogorov complexity and high 
generalization capability. Neural Networks, 10(5):857-873, 1997. 



J. Schmidhuber. Algorithmic theories of everything. Report IDSIA-20-00, |quant- 
ph/00lll22|, IDSIA, Manno (Lugano), Switzerland, 2000. 



J. Schmidhuber, J. Zhao, and M. Wiering. Shifting inductive bias with success- 
story algorithm, adaptive Levin search, and incremental self- improvement. Machine 
Learning, 28:105-130, 1997. 

C. E. Shannon. A mathematical theory of communication. Bell System Technical 
Journal, 27:379-423, 623-656, 1948. Shannon-Fano codes. 

R. J. Solomonoff. A formal theory of inductive inference: Part 1 and 2. Inform. 
Control, 7:1-22, 224-254, 1964. 

R. J. Solomonoff. Complexity-based induction systems: comparisons and convergence 
theorems. IEEE Trans. Inform. Theory, IT-24:422-432, 1978. 

R. J. Solomonoff. Applications of algorithmic probability to artificial intelligence. 

In Uncertainty in Artificial Intelligence, pages 473-491. Elsevier Science Publishers, 
1986. 

V. Strassen. Gaussian elimination is not optimal. Numerische Mathematik, 13:354- 
356, 1969. 

P. M. B. Vitanyi and M. Li. Minimum description length induction, Bayesianism, and 
Kolmogorov complexity. IEEE Transactions on Information Theory, 46(2):446-464, 
2000. 

A. K. Zvonkin and L. A. Levin. The complexity of finite objects and the development 
of the concepts of information and randomness by means of the theory of algorithms. 
RMS: Russian Mathematical Surveys, 25(6):83-124, 1970. 



